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Self-interacting dark matter resolves the issue of cuspy profiles that appear in non-interacting 
cold dark matter simluations; it may additionally resolve the so-called "too big to fail" problem 
in structure formation. Asymmetric dark matter provides a natural explanation of the comparable 
densities of baryonic matter and dark matter. In this paper, we discuss unique indirect detection 
signals produced by a minimal model of self-interacting asymmetric scalar dark matter. Through 
the formation of dark matter bound states, a dark force mediator particle may be emitted; the decay 
of this particle may produce an observable signal. We estimate the produced signal and explicitly 
demonstrate parameters for which the signal exceeds current observations. 



I. INTRODUCTION 

A number of scenarios have been put forth for ex- 
plaining why the amounts of dark matter and ordinary 
matter are relatively close to each other, within one 
order of magnitude. One of the popular approaches 
is to consider dark matter with a conserved particle 
number and a particle-antiparticle asymmetry related to 
matter-antimatter asymmetry [THU]. (For a recent re- 
view, see, e.g., Ref. [Hj.) Such asymmetric dark matter 
does not annihilate at present, and, as long as it is sta- 
ble, no indirect detection signals are expected in gamma 
rays or neutrinos. At the same time, some inconsisten- 
cies between numerical simulations of cold dark matter 
(CDM) and the observations hint at the possibility of self- 
interacting dark matter Indeed, interactions of 
dark-matter particles can facilitate the momentum trans- 
fer and angular momentum transfer in halos, hence cre- 
ating cored rather than cuspy density profiles in both 
dwarf spheroidal galaxies and in larger halos. There is a 
variety of particle-physics candidates for self-interacting 
dark matter, which include, e.g., hidden-sector particles 
with gauge [30], or Yukawa interactions [S^, as well as 
non-topological solitons with a large enough geometrical 
size 

In this paper we will show that, if dark matter is 
both asymmetric and self-interacting, then its detection 
in gamma rays is possible. Although asymmetric dark 
matter particles do not annihilate, their interactions can 
result in emission of quanta of the field that mediates 
the self-interaction. We will investigate several methods 
of producing these quanta and decay paths; on the pro- 
duction side, we particularly consider emissions occuring 
in elastic scattering (as bremsstrahlung) , or in the events 
of two interacting particles forming a bound state. The 
latter is plausible because a number of models have em- 
ployed attractive self-interaction, such as, for example 
Yukawa fields. 

To illustrate the possibilities of indirect detection, we 
will consider a fairly generic model of scalar dark matter 



S interacting by means of exchange of some lighter scalar 
field (7, both of which are singlets of the Standard Model 
gauge group. The mediator field cr can have a nonzero 
mixing with the Higgs boson, and, therefore the a boson 
can decay into photons and other Standard Model par- 
ticles, even if its coupling to Standard Model particles 
are otherwise highly suppressed. This decay ultimately 
produces the signal detectable by gamma-ray telescopes. 

The paper is organized as follows. In section|ITj we dis- 
cuss the relevant features of the particle physics model 
under consideration. In section |III[ we consider the 
bound states formation and the prospects of indirect de- 
tection of dark matter forming bound states. Finally, 
in section |IV[ we discuss possible signals from bremm- 
strahlung of self-interacting dark matter. 



II. THE MODEL 

We begin this section by introducing the minimal par- 
ticle physics model that we will use. Following this, we 
discuss properties of the dark matter halo within the 
Milky Way. We then discuss the relevant constraints on 
the parameters present in the model in general, although 
detailed discussion of the implementation of these con- 
straints is contained in later sections where the signal 
from bremsstrahlung and bound state formation is ex- 
plicitly calculated. Then finally we discuss the decays of 
the dark force mediator. 



A. Dark Sector 

We begin with a specific model of the relevant parti- 
cle physics. We supplement the Standard Model with a 
complex scalar SUc{^) x SUl{'^) x ^7^(1) singlet S; we 
also introduce a global ?7s(l) symmetry under which 

S^e^'S S^-^e-"^S'^. (1) 
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Without a loss of generality we may assume S particles 
carry unit J7s(l) charge. Due to charge conservation, the 
S particles are completely stable. 

We will assume that dark matter is composed of the 
S particles, and the correct abundance is generated in 
some process similar to or combined with baryogenesis. 



We do not assume that dark matter is necessarily a ther- 
mal relic. To make this dark matter self-interacting, we 
introduce an additional scalar field a which is a singlet 
under all the gauge symmetries, as well as C/s(l). The 
most general potential, after the Standard Model gauge 
symmetry is spontaneously broken, is 



V = M'h^a + mlS^S + -^a^ + ^{h'>f + + A,n^ + A^gS^Scr + A,h{h'fa + A^.a^h^ + A^sS^ Sh'' 



— a 
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(2) 



The Higgs field and the a field carry identical quan- 
tum numbers and therefore mix. The true mass eigen- 
values are 



m 



1:2 
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and the eigenstates are 



= cos{0M/2)h° + sin(6'Af/2)CT 
02 = - sm{eM/2)h° + cos(6'm/2)(t 



where the mixing angle is 



tanf 



M 
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(4) 
(5) 

(6) 



We will require that the mixing between the Higgs field 
and the a field be small; this can be accomplished by 
setting the free parameter M appropriately. Then we 
may speak of the a fields and Higgs fields as approximate 
mass eigenstates with masses irih and ma respectively; 
this allows the mass nia to be small even though no light 
scalar boson has been observed. 

We assume that any interactions between the dark sec- 
tor particles {S, , and tr) and the particles of the Stan- 
dard Model are highly suppressed. 

Finally, we note that the relevant unitless coupling to 
describe the Yukawa interaction between the S and a 
particles is a = A^g/l&nrrvg. This can be estabhshed in 
two ways. Regarding bound states, it is well-known that 
the Bethe-Salpeter equation reproduces the ground state 
energy of the hydrogen atom. Therefore, one may deter- 
mine a by setting the binding energy of the lowest bound 
state, ^^g/10247r^m|, equal to a^ms/4, which includes 
the correction for identical particles. Secondly, one may 
consider the non-relativistic limit of two particle scat- 
tering. We recall that quantum-field-theoretic wavefunc- 
tions include a normalization factor of \/y/2ms for each 
S particle. Therefore, the relevant prefactor before the 
overlap integral for one particle exchange is 47ra = ^^5/ 
4m|., which again gives a = A^g/WiTm^. If one instead 
defines a = A^g/ AiTm^g, as in 38] for example, then addi- 
tional factors of 4 must be introduced in other equations, 
e.g., the bound state mass. 



B. Dark Matter in the Milky Way 

Let us now discuss the assumptions that we will make 
regarding the properties of dark matter in the Milky Way 
halo. First, we assume that the correct abundance of S 
particles is determined by some process that is similar 
to baryogenesis or related to baryogenesis, as in models 
reviewed in Ref. [5T]. The absence of antiparticles in 
today's universe eliminates the possibility of a signal from 
SS'^ annihilation. 

We use the Navarro-Frenk- White profile 39J to approx- 
imate the spatial mass distribution of dark matter. 



Po 

{r/R,){l + r/R, 



(7) 



We do not expect this profile to be accurate near the 
center of the galaxy; indeed, one of the motivations of 
self-interacting dark matter is to remove the cusp present 
at r = in the NFW profile. Therefore, we will cut off 
our intergrals at scales of 1 kpc. We emphasize that our 
results are not dependent on the sharp cusp present in 
the NFW profile. The parameters po and Rs are related 
to the virial mass, virial radius, and concentration by 



Rg 
Po 



I vir 
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ln(l + C)-C/(l + C) 47ri?3' 



(8) 



For the Milky Way, we use the parameters Mvir = 1.0 • 
IQi^Afo, rvir = 258 kpc, and C = 12 gD]. This gives 
Rs = 3.4 • 10^6 GeV"^ and pn = 1-4 • 10-42 GeV*. 

In calculating the cross sections for bremsstrahlung 
emission of a particles and bound state formation, we 
will need to average over the relative velocities of the 
particles. Therefore, we need the velocity distribution 
P{v{r)) as a function of the distance from the center of 
the galaxy. If the dark matter has virialized, then its av- 
erage circular velocity should decrease near the galactic 
center, except for a small region near the supermassive 
black hole. However, the dark matter radial velocity pro- 
file and dispersion are currently unknown. 
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Because of these uncertainties, we will instead use a 
Maxwellian distribution with the effective temperature 
Tcff chosen such that the average velocity is 220 km/s. 
We note that simulations support the assumption of a 
locally Gaussian velocity distribution even for cold dark 
matter |41) . and the isothermal approximation is better 
for self-interacting dark matter [13]. The velocity 
distribution for two non-relativistic S particles is 

P{vi,V2) dvi dV2 

V^TTi off / 

(9) 

In terms of the total velocity = fi + V2 and the 
relative velocity v-^^i = 'I'l ^ i^2j the distribution is 

(47r)2 f ms \^ 
P{vrci, vt) dvrci dvT = — ^ I 2j^Ts ) 

We integrate over the total velocity to find the relative 
velocity distribution (in a reference frame at rest with 
respect to the Milky Way). 

... 

Because the 5 particles are moving non-relativistically, 
this distribution also applies to their center of momen- 
tum frame. We observe that this is peaked at slightly 
larger velocities than the velocity distribution of a single 
particle. 

C. A General Discussion of Constraints 

Thus far, we have introduced a model which provides 
a viable dark matter candidate. We introduced several 
parameters in the Lagrangian describing our model (e.g., 
the self-interaction coupling). These parameters cannot 
be set arbitrarily; there are numerous constraints they 
must satisfy, from both astrophysics and particle physics. 
In this subsection, we will give only a general discussion 
of these constraints; the specifics of how the constraints 
are implemented will be discussed when particular values 
for the coupling constants are chosen, which will be done 
seperately for bremsstrahlung emission and bound state 
formation. Since we aim to demonstrate that this model 
produces an observable indirect detection signal, we de- 
mand that it satisfy all experimental constraints except 
those from indirect detection experiments. 

First, we require that this model make only insignifi- 
cant modifications to the branching ratio for the decays 
of the Higgs boson. We forbid the decay -> SS^ by 
requiring ms > rrii^o /2 « 63 GeV, using the recent Higgs 
mass measurements [Hj, [JS]. The decay — )■ aa can 



be arbitrarily suppressed by taking Aa-h to be sufficiently 
small; this parameter is not used elsewhere in our analy- 
sis. We also demand that the mixing angle 9m be small 
enough that the apparant branching ratio for ^ a is 
less than the branching ratio for the h'^ — > 77 decay. 

There are many well-known bounds on the self- 
interaction cross section of dark matter. As explained 
in j37J, these constraints more appropriately restrict ax, 
the momentum transfer cross section. (For identical par- 
ticles, the closely-related viscosity cross section should be 
used instead |46j . In the limit msv/m^ 3> 1, which will 
be valid for our parameters, these differ by 0(1) 
and so we will ignore this complication.) The bullet 
cluster bound requires ass/ms ^ .7 cm^/g [37], and 
bounds from the evaporation of galactic halos favor crss/ 
™s ^3 -1 cm^/g |35j. These bounds appear to be in con- 
flict with the prefered range to eliminate cuspy profiles, 
.56 cmVg < ass/ms < 5.6 cmVg [21123 . 

However, because these bounds affect vastly different 
scales, they may be resolved by considering a velocity 
dependent into the cross section, as naturally arises in 
the Yukawa exchange of a light boson [JS], [37]. Fur- 
thermore, such a cross section may additionally solve the 
"too big to fail" problem [34l Ell gg. For an attrac- 
tive Yukawa potential, as we have introduced above, the 
bounds are consistent for Yukawa interactions provided 
that the masses ms and m„ satisfy particular relations 
given in [50], [51]. The precise constraint is a function of 
Wmax = \/'2ama-/TTms, the velocity at which var peaks 
at a transfer cross section equal to cr™"'^ = 22.7/m^. 

Additional bounds on the self-interaction cross sec- 
tion arise from observations of halo ellipticity, as intro- 
duced in [S2], although these bounds are quite model- 
dependent. Yukawa couplings are discussed in [3T], which 
uses the observed elliptical shape of the dark matter halo 
of galaxy NGC 720. If the self-interaction between the S 
particles is too strong, the energy transfer from these col- 
lisions makes the halo spherical instead of elliptical. Ref- 
erence [31 presented analyses for masses up to 4 TeV; 
however, we will consider masses above this. Further- 
more, as has been noted by [3^, these bounds may in 
fact be somewhat weaker due corrections from the triax- 
ial distribution of dark matter outside of the core; how- 
ever, as they note, more detailed simulations are required 
to firmly establish this conclusion. We discuss these is- 
sues in more detail in appendix |A| in which we extend 
the halo ellipticity bounds to the relevant mass range. 

Direct detection experiments such as XENONIOO [53] 
and CDMS 'ST have set an upper bound on the cross 
section for the interaction between S particles and nucle- 
ons; because this interaction occurs through the exchange 
of a Higgs boson, this constrains Ahs- The stability of 
neutron stars generally imposes stronger constraints on 
AhS [55H55] but these constraints do not apply to scalar 
dark matter with masses at the TeV scale or above [59| . 
We will not use A^s in our analysis; therefore, it can be 
set arbitrarily small. These constraints can also constrain 
the quartic interaction between self-interacting dark mat- 
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ter [BU]; we may also set this aritrarily small because it 
will not be used in our analysis. We note that while we 
can arbitrarily suppress the S-nucleon interaction which 
occurs through the exchange of a Higgs boson, there is 
an additional diagram in which the S boson emits a a 
boson, which turns into a Higgs boson via mixing, which 
is then absorbed by the nucleon. Although we are not 
free to arbitrarily suppress this diagram, as one might 
expect, this cross section is beneath current direct de- 
tection limits; we discuss it in more detail in Appendix 

m 



D. Decays of the Dark Force Carrier Particle 

From the constraints discussed above, we have seen 
that the mass mo- must be relatively small. However, 
these dark force mediator particles are not necessary sta- 
ble, and their decays can potentially produce detectable 
signals. Later in the paper we will discuss how these a 
bosons are produced (for example, though bound state 
formation or bremsstrahlung) ; in this section, we will 
simply discuss their decays irrespective of their produc- 
tion. Due to the nonzero a-Higgs mixing, a a particle 
has the same decay modes as the Higgs boson, provided 
that they are kinematically allowed. The amplitudes are 
suppressed by the cr-Higgs mixing parameters. Since the 
mass Too- must be small, we consider the decays a — ^ 77 
and fj — e+e~. 

For Too- ~ a few MeV, the dominant decay is u — > 



The decay rate in the rest frame of the cr boson is 



gl^mlma sin^(0A//2) 



327rTO^ 



1 - 



Ami 



3/2 



(12) 



where gw is the weak coupling constant. 

If m„ < 2TOe, the decay a — > e'^e~ is kinematically 
forbidden, and instead the dominant decay is cr —>■ 77. 
In the cr particle's rest frame, the decay rate is 



sm 



9 9 
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We note that in a reference frame in which the cr 
particle is moving with speed v, its gamma factor is 
7 = (1 — w^)^/^ and its hfetime is t = 7/r. 

For the product of decay to be observed, the parti- 
cles must decay in flight before they travel the distance 
« 8 kpc separating Earth from the galactic center. For 
our choice of parameters, the mean distance travelled in 
the Milky Way's rest frame before decaying is signifi- 
cantly shorter than this distance. We also note that we 
do not need to take into account scattering when calcu- 
lating decays of cr bosons because the collisions are rare 
(see Appendix [C| . 

To determine the signal produced by our model, we 
also need to know the width of the energy distribution 
of the decay products. The above decays are two body 
decays; therefore, in the rest frame of the cr particle the 
energy spectrum of the decay products is a sharp line at 
TO0-/2. This energy spectrum must be boosted into the 
Milky Way reference frame, in which the cr particles are 
moving with speed v = ^JE"^ — rn^/Ea- Because the a 
boson is spinless, the energy distribution is flat. 

For decays to an electron and positron, the energy dis- 
tribution is 



^{El-ml){l-Aml/ml 



(16) 



for Ep between the values of 



_E„ ^ ^{El-ml){l- Ami/ml) 

Similarly, for decays to two photons, the energy distri- 
bution is 



P(E,) 



1 



^El - ml 
for Ep between the values of 



Ea J El - ml 



III. DARK MATTER BOUND STATES 



(17) 



In this equation, a w 1/137, gw is again the weak 
coupling constant, Nd is the number of color states of 
the particle in the loop, and is this particle's electric 
charge. The dominant contributions to the loop will be 
from electrons, up quarks, and down quarks, for which 



^^ = -r(l + (l-r)/(r)) 
where r = Am'f/mf^, and 




T > 1 
T < 1 



(14) 



(15) 



It has previously been observed that many models 
of self-interacting dark matter, including supersymmet- 
ric models, permit the existence of dark matter bound 
states [38]. The same reference notes that the decay of 
emitted force carrier particles could, in theory, produce 
a signal for indirect detection experiments. Therefore, 
we will begin by explicitly calculating the produced sig- 
nal for the above asymmetric scalar dark matter model. 
We will establish that it is indeed possible to produce a 
signal above current observational bounds, establishing 
the possibility of indirect detection of asymmetric self- 
interacting dark matter. However, we will further show 
that the limit a <C 1, as taken in [35* , does not produce 
a detectable signal. 
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A. Choice of parameters 

We begin by discussing in more detail the constraints 
that our parameters must satisfy. To facilitate the for- 
mation of bound states, we desire a large coupling Aa-s', 
we will choose a = A^g/167rm| = 2. Although calcu- 
lations in the strongly-interacting regime are notorously 
difficult, the astrophysical bounds determined in [3T] and 
[Sn] may be extrapolated to these regimes; in both refer- 
ences, the transfer cross section used includes corrections 
for the strongly-interacting regime, as the authors note. 
We will also consider a — I and show that this is not 
sufficient to produce an observable signal. 

First, we ensure that our parameters are consistent 
with ellpitical halos; as discussed in 31j, the restriction 
on a becomes weaker as mg- is increased. In Appendix 
[A] we have determined the minimum TOo- for which we 
may consistently choose a = 2 as a function of mg. We 
choose ms = 4 TeV, which requires that we choose ^ 
30 MeV; we will take = 40 MeV. Similarly, for 
ms = 4 TeV and a = 1, we must satisfy rrio- }^ 20 MeV; 
we choose rria- = 25 MeV. 

We must also ensure that our parameters are consis- 
tent with the astrophysical data. The velocity for which 
vax = crip"-^ is Wmax = \/2amcr/7rm5; this is 1100 km/s 
and 2000 km/s respectively for the two sets of parame- 
ters above. To be consistent with astrophysical data, we 
then must have ct^'^^/tos < 100 GeV"^ 50J; the above 
numbers correspond to .2 GeV~^ and 9 GeV~^. 

In order to produce a bound state, a real a particle 
must be emitted; therefore, we must also have ma- <^ B, 
where the binding energy B = a'^ms/4:. For the first set 
of parameters chosen, the binding energy is 4 TeV, and 
for the second set of parameters, B = 1 TeV. For both, 
rua < B. 



B. Production of Bound States 

The rate of formation of bound states, neglecting 
charge depletion, is given by 



dNBS 
dt 



(18) 



where ns(r) — p{r)/ms is the number density of S parti- 
cles and ctbs is the cross secton for bound state formation. 
Because the S particles do not escape to infinity, this can- 
not be approximated using the Born approximation. In 
Appendix [P] we present a calculation of this cross section 
as a function of the relative momentum |p| = /xurci where 
fi is the reduced mass of the system. (This cross section 
is calculated by adapting the derivation for positronium 
formation given in [61] to the case for the exchange of 
a spinless boson.) The distribution of relative momen- 
tum of the incoming particles can be found from equa- 
tion (11) and is given by equation (D19). Averaging the 



above equation over the relative momentum gives 
dN^s 



dt 



ns{rfdV 

[[^^a{\p^J)P{\p\)d\p^J. (19) 



The first set of parameters discussed above corresponds 
to a cross section of 4.11 • 10~^ GeV, which gives the rate 
dNBs/dt = 2.1 • 10" GeV. In one year, 9.8 • 10^^ bound 
states are formed, which means that during the lifetime 
of the Milky Way, 1.3 • 10^^ would have formed. This is 
indeed negligible in comparison to the total number of 
S particles between 1 kpc and 8 kpc, which is 7.2 • lO^'^. 
This justifies our neglect of charge depletion. 

If we decrease a to 1, then the cross section drops by 
two orders of mag nitude, to 5.76 • lO""* GeV"^ The 
rate is also two orders of magnitude smaller, dN^s /dt = 
2.9 • 10^^ GeV. Again, we may neglect charge depletion. 



C. a Boson Production and Decay 

The S particles do not interact electromagnetically. 
Therefore, when a bound state is formed, the excess en- 
ergy is carried off by a light a particle. Although the 
binding energy is large enough that a Higgs boson could 
be emitted instead, the a particle dominates because it 
is lighter and has a stronger coupling to the S particles. 
Due to the nonzero tJ-Higgs mixing, this a particle has 
a non-zero probability to transform into a Higgs boson 
which then decays. Given our choices for m^, the domi- 
nant decay is cr —> e+e^. 

In the rest frame of the Milky Way, the a particles 
will have a typical energy equal to the binding energy; 
the additional energy the a particle may carry from the 
kinetic energy of the non-relativistic S particles is negligi- 
ble. Using equation ( 12 1, we find that the lifetime of the a 
boson, in the Milky Way's rest frame, is 1.46 • 10^^ GeV, 
or 9.61 s, for the first set of parameters. The distance 
that they travel before decaying is 10^ m, which is sig- 
nificantly less than both the distance from the galactic 
center to the solar system and the mean free path calcu- 
lated in Appendix [C] For the second set of parameters, 
the lifetime is 9.36 • 10^^ GeV = 6.17 s. 



As explained in section II D the resulting electrons and 
positrons have a flat energy distribution; their spectrum 
is 



dN, 



1 



BS 



dE^ dt ^(B2-m2)(l-4m2/m2) 



dt 



(20) 



These electrons and photons have typical energies on 
or just below the TeV scale; they lose energy through 
synchrotron radiation and inverse Compton scattering 
rapidly, within about 1 kpc Therefore, few of these 
particles will be observed near Earth. 
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D. Scattering of High Energy Electrons and 
Positrons 



There are three sources of background photons: 
the cosmic microwave background (CMB) radiation, 
starhght, and the starhght reprocessed by dust (including 
the extragalactic background hght, which is the starhght 
re-emitted by dust outside Milky Way) . Outside the cen- 
tral molecular zone, the cosmic microwave background 
radiation dominates the photon number density [63J. If 
the cross section exhibited a rapid growth with energy, 
the higher-energy optical photons could have played a 
role, but this is not the case for the IC cross section. 
For the signal from distances between 1 kpc and 8 kpc 
from the galactic center, one may safely neglect scatter- 
ing from photons other than CMB photons. Due to the 
smallness of the IC mean free path, other propagation 
effects are not significant. 

In Appendix [E] we calculate the final energy distri- 
bution for a scattering photon using the Klein-Nishina 
cross section; we then average over the appropriate en- 
ergy distribution for the electrons and positrons produced 
by dark force mediator decays. The rate of the produc- 
tion of these photons is the rate of production of the 
high energy fermions themselves. The end result of this 
calculation is described by equation (E15l, which gives 



the number of photons produced by this process per unit 
energy per unit time; from this, we can find the flux of 
gamma rays at the solar system. 

The production of dark force mediator particles results 
in an isotropic flux of these particles about the galactic 
center; similarly, we expect the flux of their decay prod- 
ucts and the scattered photons to be isotropic about the 
galactic center. Therefore, the photon flux per unit area 
is well-approximated by an equivalent point source at 
the galactic center. We can find the average flux per 
unit area, per unit solid angle by further dividing by 27r, 
since the signal will appear to come from the hemisphere 
centered on the galactic center. We note that this is an 
average; as a function of solid angle, we expect the sig- 
nal to be greater near the galactic center and less further 
away from it. We also note that this is only an approx- 
imation to the true diffuse flux, meant to demonstrate 
that a detectable signal is possible. (As discussed above, 
we neglect the contribution from the galactic center itself 
and include only the contributions from decays outside 
the inner kiloparsec.) If such a signal were to be observed, 
more careful analysis should be done before attempting 
to fit this scenario to the data. 

Furthermore, we note that the production of the dark 
force mediator bosons a scales as the density squared; 
this increases as one approaches the galactic center. Since 
the sigma bosons only travel 10^ m before decaying into 
the fermions which scatter the CMB photons, the signal 
is dominated by the innermost region we consider. Since 
we have cutoff our calculation at an inner radius of 1 
kpc to avoid the known cusp in the NFW profile, the 
signal comes predominantly from the region near this cut. 




0.5 1 5 10 50 100 500 

Energy of Scattered Photons E (GeV) 



FIG. 1. Depending on the model parameters, the signal can 
range from undetectable to already excluded. The signal is 
shown for two values of a, as shown in the legend, and for 
ms = 4 TeV, = 40 MeV for a = 2, m„ = 25 MeV 
for Of = 1. For comparison, the data from Fermi LAT space 
telescope are also shown |64) . 



Therefore, the point source approximation is better than 
one may naively expect. (We remind our reader that due 
to this cutoff, this calculation produces an approximate 
lower bound on the signal strength.) 

We find that the average flux over the hemisphere cen- 
tered on the galactic center, neglecting the galactic center 
itself, is 



^ ^ diV-y,tot ^ 1 

dE dt ' 27r St ' 47r(8 kpc)^ ' 



(21) 



To compare with the the sensitivity of the Fcrmi-LAT 
Gamma Ray Telescope, we evaluate i?^$; this function 
is plotted in Fig. [l] The signal for a = 1 is peaked 
at a lower energy and falls off more sharply, as we would 
expect because the binding energy is smaller. We see that 
a = 2 produces a signal that is one order of magnitude 
larger than the values measured by Fermi-LAT, but the 
signal produced by a = 1 is two orders of magnitude 
too small. Therefore, we conclude that sufficiently large 
couplings may produce a detectable signal. This suggests 
that WIMPonium models [3S], which assume a <C 1, 
will not produce a detectable signal through bound state 
formation. 

One can also show that the resulting signal is rather 
insensitive to the precise value of m„, provided that a — 
e+e~ remains the dominant decay. However, this param- 
eter is highly constrained by the astrophysical bounds 
discussed in section [II C[ Finally, we also note that this 
signal depends relatively weakly on the cutoff we imposed 
to avoid the center cusp of the NFW profile. If we cut off 
the integral at 1 pc instead of 1 kpc, the signal would only 
be about 20 percent greater, although the point source 
approximation would be more accurate. 



7 



E. Possibility of a Positron or Electron Excess 

We will briefly discuss the possibility of producing a 
detectable positron or electron signal within this model. 
This is a particularly interesting question because of the 
positron excess observed by PAMELA [65i, which was 
confirmed by Fermi-LAT [55] and more recently AMS- 
II 67]. In order to travel from the galactic center to the 
solar system relatively unimpeded, the fermions would 
need to be lower energy that those discussed above, which 
lost significant energy due to inverse Compton scattering. 
Energy loss due to inverse Compton scattering is some- 
what suppressed for energies on the GeV scale; therefore, 
we will briefly discuss the difficulties of producing elec- 
trons and positrons on this scale. 

To produce a significant number of dark force mediator 
particles, we desire to keep a relatively large in order to 
maintain a large cross section for bound state formation. 
However, the energy of the fermions produced by the 
decay of the dark force mediator depends only on the 
mass of the S particles and a. Therefore, to produce 
10 GeV-scale positrons and electrons, we must decrease 
ms to be on the scale of a few hundred GcV. This is 
below the scale typically discussed in the WIMPonium 
literature. 

Naively, these parameters appear to run into difficul- 
ties with the halo ellipticity bounds such as in ^T\; for 
example, ms = 100 GeV with a = 2 appears to re- 
quire = 232 MeV, for which the dominant decay is 
to muons instead of e"'"e~. (However, it should be noted 
that the analytic approximation for the cross section be- 
gins to break down at ^ 100 MeV.) This appears 
to eliminate the possibility of an observable electron or 
positron excess. 

However, a more detailed analysis of the bounds by 
[35] suggests that these bounds should be about an or- 
der of magnitude weaker. In Appendix [^ we have 
parametrized this uncertainty with the parameter F, 
which is one if the considerations of p35] are ignored. If 
one assumes F ^ .1, then a small region of parameter 
space remains which is consistent with the halo elliptic- 
ity bounds and is small enough (50 to 80 MeV) that 
the decay to e+e~ dominates. 

The analysis proceeds as above, up to the point where 
one calculates the inverse Compton scattering. For these 
lower energy electrons and positron, we do not expect in- 
verse Compton scattering to be a significant effect. How- 
ever, other effects can influence the shape of the spec- 
trum observed at Earth; for example, we must be partic- 
ularly concerned with positron annihilation, which will 
generically decrease the detected positron fraction. We 
do expect these particles to lose energy due to bremm- 
strahlung. A detailed analysis could be run using cosmic 
ray propagation software such as GALPROP. 

However, we believe that it is unlikely that the re- 
sulting spectrum could be tuned to reproduce the ob- 
served positron excess observed in PAMELA [55] , Fermi- 
LAT [55], and AMS-II [57]. The energy spectrum of 



fermions produced through the decay of dark force me- 
diator particles is flat, and while this spectrum will no 
doubt be modifled by a detailed analysis of propagation 
from the galactic center to the solar system, we think the 
resulting E^^ is unlikely to be as flat as that observed by 
PAMELA. Furthermore, the positron excess extends to 
higher energies beyond that which can be accomodated 
by our model. Hence, we conclude that such models are 
unable to account for this observed excess. 

Finally, we note that it may be possible to adjust the 
parameters so that an excess of positrons or electrons 
above PAMELA'S observations is produced, although 
again careful analysis of the propagation of said fermions 
would be necessary. If such an excess can be produced, 
we would expect it decrease at or before the TeV scale, at 
which point the spectrum would be limited due to inverse 
Compton scattering energy losses. Since no such behav- 
ior is observed in PAMELA'S spectrum, one would trans- 
late this into bounds on the parameters. However, since 
the parameter space in which such a signal is potential 
possible is already quite small, any resulting constraints 
(if any) would be quite weak. 



IV. BREMSSTRAHLUNG EMISSION OF DARK 
FORCE MEDIATORS 

In this section, we will discuss the signal produced by 
bremmstrahlung emission of a cr boson which decays to 
photons. First, we choose our free parameters consistent 
the constraints discussed in III CI Then we calculate the 
rate at which a bosons are produced via bremsstrahlung 
near the galactic center. Finally, we determine the flux 
at the Earth and compare with observations by the IN- 
TEGRAL experiment. We will see that while the flux 
exceeds the flux produced via bound state production, 
the resulting signal is significantly beneath observational 
limits, due to the substantial background at lower ener- 
gies. 



A. Choice of parameters 

In bremsstrahlung emission, the emitted a boson will 
carry an energy comparable to the kinetic energy of the 5' 
particles. Since these have a velocity of order 10~^, this 
means that the typical energy scale of bremsstrahlung 
emission will be six orders of magnitude below mg. We 
will require mg. ^ v'^ms/2 because we do not want an 
additional suppression from the difficulty of emitting real 
a bosons. This will infiuence the implementation of the 
constraints discussed in section II C (We note, however. 



that due to the contribution of the tail of the relative 
velocity distribution, we do not necessarily expect a sharp 
cutoff at the average kinetic energy.) 

The first constraint comes from Ref. [50], which de- 
termines the condition for astrophysical observations to 
be consistent with dark matter whose self-interaction is 
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described by a Yukawa pot ential. The p recise constraint 
is a function of Wmax = \/'^OLm„l'Km,s, the velocity at 
which va^p peaks at a transfer cross section equal to 
(j-rnax _ 22.7/m^. If Umax ^ 10 km/s, then the astrophys- 
ical constraints are consistent if 22.7 /m?^ms ^ 35 cip? / 
g = 16000 GeV^'^ 5Tj; we will verify that we satisfy this 
condition below. Combining this with <C v^ms/2 
reveals that we must satisfy 



22.7-4 



■ 16000 GeV" 



1/3 



1.3 TeV. 



(22) 



Let us choose ms — 10 TeV and mg. — .5 MeV. 

Next let us discuss the bound from halo ellipticity; this 
will constrain the coupling A^ys- We extended the anal- 
ysis of Ref. [3T] to ms — 10 TeV in Appendix [Aj this 
showed that for = .5 MeV and ms = 10 TeV, we 
prefer to take a — A^g/Wirm's ^ .93 although this 
may be loosened somewhat. If we choose to saturate 
this bound, we find A^s = 68 TeV. These values give 
Vmax ~ 50 km/s, self-consistent with our initial assump- 
tion that Umax ~ 10 km/s. 

The bremsstrahlunged a boson can be emitted by ei- 
ther of the S particles; however, it can also be emitted by 
the a boson exchanged between the S particles. These 
diagrams involve the coupling Ac, which is thus far un- 
constrained. In order to enhance the signal, we will sat- 
urate the perturbativity bound, taking A^- = 3.5 MeV. 
On the other hand, we will also consider A^ — 0, which 
is equivalent to neglecting the two diagrams on the right 
of Fig. |2] While this certainly won't help to increase our 
signal, the properties of the signal will be qualitatively 
different in the two cases in interesting ways. 



B. Production of a Bosons through 
Bremsstrahlung 

Next, we must know the cross section for 
bremsstrahlung emission of a soft a boson, which 
involves evaluating the 10 diagrams shown in Fig. [2] 
The derivation of this cross section, including averaging 
over the relative velocity of the incoming particles, is 
contained in Appendix [Fj For the parameters given 
above, the cross section is a = .0108 GeV^^. This 
is the same order of magnitude as the a = 2 cross 
section for bound state formation; we note, however, 
that we did not have to increase the coupling a into the 
non-perturbative regime in order to reach this value. In 
general, as we would expect, the bremsstrahlung cross 
sections are indeed large in comparison to the bound 
state formation cross section. 

We shoud note that this cross section does not include 
any enhancement due to Sommerfeld factors; this contri- 
bution will be discussed later. Also, although naive esti- 
mates would suggest a large enhancement, in this regime 
the Sommerfeld factor may be unreliable and a proper 
resummation suggests that any enhancement is at most 



0(1) to 0(10)155]. This is discussed in somewhat more 
detail after the calculation of the cross section in Ap- 
pendix |Fj 

The rate of production of bremsstrahlung a bosons is 



dt 



ns{r) 

^rcl'^brcm {Vrcl)P{Vrcl) dV dv^d (23) 



where ns{r) = p{r)/ms is the number density of dark 
matter S particles, and pir) is given by equation ([7|. 
We have also averaged over the relative velocity of the S 
bosons, and the integration extends from 1 kpc to 8 kpc, 
the distance from the solar system to the galactic center. 
For the given parameters, dNa/dt = 3.21 • 10^'^ GeV, 
or 4.87 • 10^"^ s^^. As might be expected, for A^, = 0, 
we find the lower rate dN„/dt = 1.45 • lO" GeV~^ = 
2.20 • 10"^^ s^^. The fact that the cross section drops by 
two orders of magnitude shows that at A^r — 2.5 MeV, 
the diagrams in which the bremsstrahlung a boson is 
emitted by the exchanged a boson dominate. Since these 
diagrams are absent for A^ = 0, we expect the signals 
produced to have qualitative differences. 

We will show that bremsstrahlung will not produce 
a detectable signal, while we found that for sufhcicntly 
large couplings bound state formation can. Since this is 
perhaps a surprising result, one may find it beneficial to 
compare with the calculation of the bound state signal 
at each step to determine why this is so. We emphasize, 
however, that such comparisons must be made carefully, 
since the bound state calculations were performed in a 
different region of parameter space. We wish to empha- 
size that for any fixed perturbative value of a, the rate of 
bremsstrahlung production will always be much greater 
than the rate of bound state formation, as one would 
expect. However, if one compares the value of dN^/dt 
found above with dN^g / dt given in the previous section, 
which are evaluated at different parameters, one finds 
that dN^s/dt is larger by about an order of magnitude, 
even though we have chosen parameters such that the 
cross sections are comparable. This is a result of taking 
ms — 10 TeV here as opposed to 4 TeV above; increasing 
ms decreases the number density ns{r). 

Next we observe that the spectrum of the emitted a 
bosons per SS — SSa event is given by 

dNcr 1 rfcTbrcm ^24^ 



dE„ 



dE„ 



As we might expect, this spectrum is sharply peaked at 
600 keV, which is on the same scale as the kinetic energy. 
The spectrum of the produced a bosons per unit time is 



dtdE„ 



ns{rfdV 



^'^brcm 7-,/ \ 7 
«rol— 77^^ P(Ui.cl) aWrcl- 



dE„ 



(25) 



C. Decay of the a Bosons and Resulting Signal 

For m„ = .5 MeV, the dominant decay mode of the 



a boson is cr — 77, which is decribed by equation ( 13 1 
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FIG. 2. These diagrams contribute to the emission of a bremsstrahlung a boson. The solid lines represent S bosons, while 
the dashed lines represent a bosons. The top line represents f-channel scattering, while the bottom line represents w-channel 
scattering. 



If we assume the mixing angle between the a boson and 
the Higgs boson is 10^'^, then the typical lifethne of the 
produced a bosons is 10^ s, during which they travel 
about 10^^ m, which is significantly less than the 10^° m 
between the galactic center and the solar system. 

The spectrum of the photons produced by the decay 
of the bremsstrahlung a bosons is given by 

where the distribution of photon energies, as a function 
of the initial a boson energies, is given by equation ( |17| . 
(The Ej dependence appears in evaluating the Heaviside 
step functions.) As we would expect, this spectrum is 
peaked around 300 keV. We note that the tail decreases 
less rapidly as A^r is decreased. As a result, the signal for 
Aa- = will be skewed torwards higher energies. 

The production of dark force mediator particles results 
in an isotropic flux of these particles about the galactic 
center; similarly, we expect the flux of their decay prod- 
ucts to be isotropic about the galactic center. Therefore, 
the photon flux per unit area is well-approximated by an 
equivalent point source at the galactic center. We can 
find the average flux per unit area, per unit solid angle 
by further dividing by 27r, since the signal will appear to 
come from the hemisphere centered on the galactic cen- 
ter. We note that this is an average; as a function of solid 
angle, we expect the signal to be greater near the galactic 
center and less further away from it. We also note that 
this is only an approximation to the true diffuse flux, 
meant to demonstrate that a detectable signal is possible 
(and we remind the reader that we are already neglicting 
the contribution from the galactic center itself). 

We find that the average flux at the solar system is 

$ = • — 27) 

Aird^ 2tt St dE^dt' ^ ' 

where d = 8 kpc is the distance from the galactic center 
to the solar system. Since we have calculated the number 



of produced tr bosons out to a radius of 8 kpc, the signal 
we calculate here comes from the hemisphere centered 
on the galactic center, which explains the 2n st. We 
note that this is the average over the hemisphere; the 
flux will be somewhat greater towards the galactic center 
and somewhat less towards the edges; however, this is a 
relatively small effect, contributing perhaps an order of 
magnitude increase as we approach the center. 

Again, with the same caveats as above, let us compare 
with the bound state case. The photon energies here are 
spread out over the scale of 100 keV, whereas the photon 
signal for the bound state production is spread over the 
scall of 100 GeV. However, a single high energy fermion 
produces about 10^ ~ 10"^ GeV-scale photons through 
scattering off of the CMB, while each a boson produced 
through bremsstrahlung produces a mere 2 photons. As 
a result, the estimated ratio of fluxes is $biom/<I'BS ^ 10^ 
or 10^. We note that since the two scenarios are in dif- 
ferent regions in parameter space, this cannot be inter- 
pretted as the ratio of actual bremsstrahlung-produced 
photons to bound state produced photons in the galaxy. 

The relevant energy scale for bremsstrahlung emission 
is on the scale of hundreds of keV, while the relevant 
energy scale for bound state emission is on the scale of 
a hundred GeV. Astrophysical backgrounds are signifi- 
cantly larger at this smaller scale; the SPI on the INTE- 
GRAL experiment records i?^<I> on the order of 1 to 10 
keV/cm^ s st for energies 20 keV and 1000 keV gl]. The 
flux of produced photons cannot be distinguished from 
this large backaground. 

The resulting signal is shown in fig. [3] as we ex- 
pect, it is about 8 orders of magnitude smaller than the 
bound state formation signal. More importantly, the 
larger A^ = 3.5 MeV signal is about 7 orders of mag- 
nitude beneath INTEGRAL'S observations. We also can 
see the qualitative difference in the signal shapes alluded 
to above; this is because for A^ = 2.5 MeV, the right- 
most two diagrams of Fig. [2] dominate, whereas these 
are absent if A^ = 0. We observe that without these dia- 
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FIG. 3. The flux of gamma rays produced by bremsstrahlung 
emission of a particles and their subsequent decay for ms = 
10 TeV, m„ = .5 MeV. 

grams, the signal is significantly smaller, but it is peaked 
at higher energies. 

We have noted above that the calculated cross section 
does not include a Sommerfeld enhancement, because 
some analysis suggest that such large factors are unre- 
liable [68] . Even if we assume that the naive Sommerfeld 
factor given by 



1 — exp(— aTr/u) 

is accurate to arbitrarily large scales, this enhancement 
is not sufficient to produce a detectable signal. For the 
parameters in the range discussed, the enhancement is of 
order 10^ or 10^, which is still too small to produce the 
seven orders of magnitude amplification required for the 
signal to be detectable. 

The signal can be increased by increasing the cou- 
plings; and indeed, as discussed in Appendix [X] there 
is some uncertainty in the halo ellipticity bounds. To 
produce a detectable signal requires increasing the cou- 
pling a to ~ 10'^, well outside the perturbative regime 
and far beyond what can be made consistent with the 
halo ellipticity bounds. It is true that is unrestricted 
by astrophysical bounds, but in order to amplify the two 
diagrams it appears in to the scale of INTEGRAL'S ob- 
servations, we would need to take A„/raa ^ IQ^ , which 
is unreasonably large. 

Therefore, we conclude that bremsstrahlung emission 
of dark force mediator particles cannot produce de- 
tectable signals, although the photon flux is generally 
significantly larger than bound state production. One 
might consider the idea that even if the signal produced 
near the galactic center is not detectable, perhaps such 



processes enhance the gamma ray or x-ray emission of 
nearby dwarf galaxies sufficiently to be observable; how- 
ever, a simple estimate reveals that this is not the case. 
Even if the signal calculated above, for the Milky Way 
galaxy, was somehow shrunk into a dwarf galaxy 40 kpc 
from us which covered a 3° by 3° patch of the sky, the 
number of counts expected in an ideal 1 m^ detector is 
of order 10~^ keV~^ s~^, which is again well below the 
background emission. 



V. CONCLUSIONS 

We have considered indirect detection signals produced 
by a minimal asymmetric self-interacting dark matter. 
Due to the f/s(l) asymmetry, the typical indirect detec- 
tion signal from dark matter annihilation is absent in this 
model. However, we demonstrated that asymmetric self- 
interacting dark matter can, in fact, produce a strong 
signal from the processes accompanying the formation of 
bound states, as has been discussed in the WIMPonium 
literature. We have found that signals are possible for 
sufficiently large couplings. This effect makes possible in- 
direct detection of asymmetric self-interacting dark mat- 
ter. The spectrum of gamma rays can help distinguish 
coUisionless dark matter from self-interacting dark mat- 
ter. We have performed explicit calculations for several 
sets of parameters; showing that for a = 2, ms = 4 TeV, 
and ma — 40 MeV the signal would be detectable. How- 
ever, we have shown that this signal is detectable only 
in the strongly interacting regime, by showing that if a 
is decreased to 1 (keeping ms constant), the resulting 
signal is not detectable. 

Then we have discussed, albeit briefly, the possibil- 
ity that this model could, in a narrow region of pa- 
rameter space, produce a detectable excess in electrons 
and/or positrons. Additionally, we have also considered 
the signal produced by the bremsstrahlung emission of 
the a boson. This was calculated for two points in pa- 
rameter space {ms = 10, mg. — .5 MeV, a — .93 with 
Aa =3.5 MeV and A^ = 0) to demonstrate two limits of 
spectrum shape. However, we have shown that although 
the flux of gamma rays can be rather large, the result- 
ing signal is actually quite small and significantly below 
backgrounds. 

This work was supported by DOE Grant DE-FG03- 
91ER40662 and by the World Premier International Re- 
search Center Initiative (WPI Initiative), MEXT, Japan. 



Appendix A: Extension of Bounds Prom Halo 
Ellipticity 

As was noted in section |II C[ one constraint on the 
self-interaction of dark matter arises from the observed 
ellipticity of dark matter halos. In this appendix, we will 
extend the result of [3T] to higher dark matter masses. As 
has been noted by |35j , these bounds may in fact be some- 
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FIG. 4. A plot of the critical coupling a = ^^5/16717715 as 
a function of 7715. Couplings below the critical coupling are 
consistent with the elliptical shape of dark matter halos. 



what weaker due corrections from the triaxial distribu- 
tion of dark matter outside of the core; however, as they 
note, more detailed simulations are required to firmly es- 
tablish this conclusion. Therefore, we will parametrize 
our uncertainty by the coefficient F\ the numerical simu- 
lations presented in [35 could be interpretted as favoring 
F ~ .1 for the particular halo model considered. 

The dark matter halo will be spherical, as opposed to 
elliptical, if collisions which change the particle velocities 
by factors of order 1 happen frequently enough. The rate 
at which these collisions occur is given by 



Tfe = y d-'v^d%^j{vx)j{v2){nsv,,xFaT)[viJvi) (Al) 

where ctt is the momentum-transfer cross section, given 
by ctt = / (Kl (dal(Kl){\ — cos(6')), and /(w) is the dark 
matter velocity distribution. The analytic fit for ctt, the 
distribution functions, and the relevant parameters for 
NGC 720 are all available in [3T]. In this reference, they 
produce plots of numerical results for mg up to 4 TeV. 
However, we will need to consider masses above this, and 
therefore, we extend their results to higher masses. We 
note that quantum corrections to the cross section be- 
come important if the limit rfigvlm^ 3> 1 is violated; 
however, all of our parameters will be in this regime. If 
not, corrections such as those discussed in [46] should be 
included. 

In particular, we extrapolate the plot of the critical 
coupling a — A^g/lGTrrrig to ms — 12 TeV. We show 
the results for rricr — .5 MeV, relevant for section IV and 
for ma — 3 MeV. We also consider F — 1 and F = .1 as 
suggested by the results of [35]. For = .5 MeV and 
ms = 10 TeV, this bound with F — 1 requires a < .93. 



However, it is substantially loosened to a < 3.6 if we 
instead take _F = .1. 

In the discussion of the signal from the formation of 
bound states, we desire a large coupling; however, we 
must remain consistent with the halo ellipticity bound. 
As discussed in [31], this bound becomes weaker as m^ 
increases. Therefore, we calculate the minimum m^ for 
which a = 2 or a = 1 is consistent with the observed halo 
ellipticity, as a function of ms ■ The results are shown in 
Fig. [5] Again we see that taking F = .1 dramatically 
weakens the bound. 



Appendix B: S-Nucleon Interaction Cross Section 



As was noted in section II C the dominant interac- 



tion between the S boson and nucleons, which is relevant 
for direct detection experiments, typically occurs through 
the exchange of a single Higgs boson. This diagram is 
proportional to A^h , which is otherwise unconstrained in 
our model, and therefore we can arbitrarily decrease this 
coupling, thus killing the signal. 

However, as we note in the text, there is another dia- 
gram which becomes dominant at sufficiently small values 
of the coupling: the 5* boson may emit a a boson, which 
transforms into a Higgs boson via mixing and couples to 
a nucleon. This diagram involves the coupling A^s and 
the mixing angle Om- We cannot take either of these pa- 
rameters to zero without eliminating the signal, although 
the mixing angle may be quite small. Thus one cannot 
arbitrarily decrease the 5-nucleon cross section; there is 
a minimum value set by this diagram. In this Appendix, 
we will show that the contribution of this diagram is in- 
deed quite small, as expected, and causes no tension with 
direct detection constraints. 

We note that the oscillation time scale, which is given 
by Tosc = 27ri?/Ar7i^, is generally many orders of magni- 
tude smaller than the interaction time scale, which can 
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be estimated by considering the overlap of the wavefunc- 
tions. Consequently, averaging over the "detector scale" 
(nucleon size), along with the source location, will sim- 
ply give a factor of 1/2. (This is in contrast with certain 
neutrino oscillation experiments, for which Tqsc may be 
large in comparison to other experimental scales, due to 
the small Am^. In our scenario, Ato^ ~ "^lO 

The S particles under consideration are generally much 
heavier than the protons; we will masses between 4 and 
10 TeV. Therefore, in the center of momentum reference 
frame the S particles will be approximately stationary, 
while the protons approach at speeds of approximately 
220 km/s. The momentum transfer is approximately 
2mpV = 1.5 MeV, which is far below the scale at which 
the nucleon form factors must be included. 

The relevant matrix element is 



where v is the vacuum expectation value of the Higgs 
boson and u, u are spinnors for the proton. This yields 



COS — — sm ' 



Because the velocities are non-relativistic, the initial en- 
ergy squared is approximately (ms -I- m„)^ « m|, which 
gives an approximate cross section 



a w ^ • — — J — — — ^ cos sm , 

167rm| to|u4 \ 2 ) \ 2 



s 




FIG. 6. The scattering of a particles on dark matter. 



Appendix C: Mean Free Path of Dark Force 
Mediator Particles 



As was noted in section [Tl D[ the mean free path for the 
a particles in the galaxy must be greater than the dis- 
tance they would travel before decaying; otherwise, con- 
stant scattering can act like a quantum Zeno experiment 
that prevents the decay. In this appendix, we present the 
calculation for the mean free path, and show that it is 
greater than the distance from the galactic center to the 
solar system for the relevant regions of parameter space. 

Since the quartic coupling Ao-s can be made aribtrarily 
small, we will assume that the scattering is dominated by 
the Sa interaction mediated by an S boson; there are two 
diagrams that contribute, which are shown in Fig. |6] 

We assume that in the lab frame, the a particle is mov- 
ing relativistically with energy Eu, while the S particle is 
moving non-relativistically with velocity v of order 10"'^. 
We do not assume any relation between E„ and the ki- 
netic energy of the S particle. Since the cross section is 
a relativistic invariant, we may evaluate it in the center 
of momentum frame, which under the above assumptions 
is attained by boosting hy P — Ecr/{Ecr + ras). Keeping 
only the largest terms, we find that the initial and final 
four-momenta in the CM frame are 



Let us consider one of the sets of parameters used 
in the bound state cross section; ms — 4 TeV and 
AaS = 20 TeV, corresponding to a = 2. For the av- 
erage effective mass of a quark, we use 3 MeV, and we 
choose 9m = 10-3. This gives a ~ 10"!^ GeV"^. All 
of our other choices for parameters give a cross section 
below this value. This is well beneath the limits from 
direct detection experiments, which are 10"'''^ cm^ or 
10"^^ GeV"^ |S3], [M]. There is a significant uncertainty 
in the contribution of the s-quark to the effective quark 
mass. Since the Higgs coupling to s is much greater than 
the couplings to u and d, even a relatively small contribu- 
tion of the sea quarks with higher masses can dominate 
the cross section. The measured s quark contribution, 
manifest as the nuclear pion-nucleon sigma term, is un- 
certain, and the resulting uncertainty in the cross section 
can be as large as an order of magnitude 70 . However, 
even at the upper edge of the range, the cross section 
does not reach the present lowest cross sections acces- 
sible in experiment. Hence, at present, direct detection 
experiments do not constrain the scenario we have con- 
sidered. 



Pa,i = (7/3™s,0,0,7/3ms) 
(7TOs,0, 0,-7/3ms) 



p'^j = jPms, j(3ms sin(6'), 0, 7/3ms cos(6')) 

Psj = il'^s, -iP-ms sin(6'), 0, -jl3ms cos{e) (CI) 

where we have used the fact that the collision is elastic. 
The matrix element is 



-iM 



ml - {ps,i - VajY ml - {ps,i + Pa,iY 
1 + cos(6') 



(Tb 



2-f'^ml V (1 + /3) (1 + /3 cos(6')) J ' 
The cross section is given by 

1 f 



(C2) 



647r2 ./ 72(1 + PYml 



1 + cos(6') 



1287r76(H-/3)4m| 7o \l + licos{e) 

2l3 + {l-li^)H{l-P)/{l + li)) 
647rTO| 76^33(1 _,_^)5 



sin(6i) de 
(C3) 
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The mean free path is ^ = {ans)~^, where ns can be 
found using equation ([t]). Since ng depends on r, the 
mean free path will also depend on r; it is the smallest 
as we approach the galactic center. Let us consider some 
typical parameters. For ms = 5 TeV, A^s = 3 TeV, 
and E„ = 1 TeV, the mean free path at 1 pc is of order 
10^^ m. If we decrease to 1 MeV, the mean free path 
increases to order 10^^ m (at 1 pc again). These values 
are all much greater than the 10^° m between the galac- 
tic center and the solar system; therefore, requiring that 
the a particles decay before reaching the solar system 
provides a stronger bound as claimed. 



Appendix D: Bound State Formation Cross Section 

In this appendix, we calculate the cross section for two 
S particles to form a bound state through the exchange 
of a bosons. We emphasize that because the S particles 
form a bound state, they do not escape to infinity, and 
therefore the Born approximation is not applicable. We 
note that, although the coupling is strong, we are in the 
classical regime, because msv/ma- ^ 1; therefore we do 
not need to include additional quantum corrections such 
as those calculated numerically in [46] . 

We will approximate the a boson as massless. The 
cross section for non-relativistic electrons and positrons 
to form a bound state through photon exchange was cal- 
culated in J61J; we adapt this derivation for scalar fields. 
The matrix element is 



■ «'z(ri,r2)dVi d^r2{2TT)S{E, ^ Ej - E,). 

(Dl) 

In this equation, ri and r2 are the locations of the two 
S particles respectively, Vf/ is the wavefunction of the 
bound state, and is the wavefunction for the two in- 
coming S particles. The factor e"*'" ''" represents the 
wavefunction of the a particle, and the sum is over the 
two 5* particles it can couple to. In this equation, the 
wavefunctions have the standard normalization in quan- 
tum field theory; however, since we are interested in the 
non-relativistic limit, let us use wavefunctions that are 
normalized to one. Then the matrix element is 



M = ~i 



^—ik-' 



— ik-r2 ~1 



2ms 

^',(ri, r2)dVi d^r2i2n)SiE, - Ef - E,) 



ri + r2 



Next we define 
R 

and write the wavefunctions as 



ri - 7-2 



*»(r2,ri) = e^'^-«vl/,(r.) 
*/(r2,ri) = e'-P-«*/(r) 



(D2) 
(D3) 

(D4) 



where Q = Pi + P2 is the total momentum of the initial 
particles. Similarly, P is the momentum of the bound 
state. After performing the d^R integral, we have 



;.(r) (^e*'-/2 + e-'^ "-/^^ ^',(r)dV 
• {2tt)'^5{E, - Ef - E„)5^{Q k P). (D5) 
The reduced matrix element is 



M 



/2 



and the differential probability is 
TV 



dW 



^6{E, 



E 



{2ttY2E„ Ami 
■ \M\'^\k\'^ d\k\dQ.d^P, 



f 



E,)S^Q^k-P) 
(D7) 



where V is the normalized volume, T is the interaction 
time, and dfl is the solid angle for the a particle. The re- 
maining integrals enforce momentum and energy conser- 
vation; we may perform them by directly imposing these 
constraints in our calculation. The transition probability 
per unit volume and unit time is 



dw 



4m| 2E„{2ttY 



If 1 



^ B then E^ ~ and this simplifies to 
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dw = 



\k\dVt 



4m| 



2(27r)2 



(D8) 



(D9) 



The differential cross section is da = dw/vrei where Vrei 
is the relative velocity of the particles in the initial state. 
We define the relative momentum by p = /iVreZ where 
/i = ms/2 is the reduced mass. |p| is also the momentum 
of one of the incoming particles in the center of momen- 
tum frame; we will now specialize to this frame. (We note 
that the cross section is Lorentz invariant, and therefore 
still applicable to other reference frames.) Then 



d<T 
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\k\dn 



Ams |p|(27r) 



(DIO) 



The free S particles do not escape to infinity; they 
exist only in the initial state. Therefore, at large r, Vl'i(r) 
must be a superposition of a plane wave and an outgoing 
spherical Coulomb wave. (Although our interaction is 
not electromagnetic, the appropriate asymptote is still a 
spherical Coulomb wave in the approximation that m„ <C 
ms-) The appropriate wavefunction to use is [7T] (also 
discussed in pT] ) 

■^^(r) = e'''^/^r(l - iC)F{iC, 1, i{pr - p ■ r))e'P"' (Dll) 

where C — A^s'ms / A\p\ms — Acts/4|p|, and F is the 
confluent hypergeometrical function. This has the same 
normalization as a plane wave. We note that the cross 
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section will be very sensitive to the ratio Aq-s/IpI ~ ^aS I 
ms as a consequence of the exponential. We adapt the 
hydrogen ground state wavefunction for ^'/(r); again, 
this is accurate in the approximation that mu is negligi- 
ble. 



(D12) 



where 77 = CIpI = this is the radius of the bound 

state. The reduced matrix element is 



F{i(, 1, i{pr ~ p ■ r)) Sr. 



(D13) 



To evaluate the integral, we differentiate the identity 



^i{p—K)-r — rir 



F{tC,l,i{pr-p-r)) 



21 -< 



All- 



+ {ri - i\p\) 
[(p-/«)2 + 772]i-< 



(D14) 



with respect to 77. The result is 

'■gz(p-^).r-„r^^j^^ 1, i{pr -p-r)) £r 



[{p - k)^ + T]^]^-< 

■^(77-i|p|)[(p-k)2 



+ [t] - i|p|)2] 



5(«,X) 



(D15) 



where x is the angle between p and k. We observe that 
g{K, TT — x) = g{~i^: x)- If the angle between k and p is 
T, the reduced matrix element is 



M = 



fcl 



Ifcl 



(D16) 



This can be evaluated numerically. The last remaining 
unknown quantity in (DlOl is |fc|, which can be found 
from the energy conservation equation 



2m + ^ ^ (2m - S) + , 

m ^ ' 2{2m-B) 



Ifcl 



(D17) 



where we have noted that in the center of momentum 
reference frame, the bound state also has momentum |fc|. 
We find the total cross section by numerically integrating 
( |D10 1. We will also average over a relative momentum 
distribution P(|p|); the total cross section is given by 



f^BS 



Als MM\ 
4ms |p|(27r) 



2P(|p|)d|p|27rsin(T)dT 



(D18) 



We note that in the non-relativistic limit the momentum 
difference of the two particles is independent of reference 
frame; therefore we can calculate -P(|p|) in any conve- 
nient frame even though we specialized to the center of 
momentum reference frame above. The total cross sec- 



tion is, of course, Lorcntz invariant. Using equation (11 1, 
we find the relative momentum distribution 



Pi\p\)d\p\ 



ms 



eff 



3/2 



-\p\^/msT,ff 



P\'d\p\. 
(D19) 



Appendix E: Scattering From CMB Photons 

When we calculated the signal produced by bound 
state formation, we found that the dark force mediator 
bosons (T decayed into TeV-scale fermions. These lose 
energy due to scattering with CMB photons, as noted 
In this appendix, we produce a calcu- 



III D 



m section 

lation of the spectrum of gamma rays produced by this 
scattering. 

First, we will show that we can neglect the energy loss 
due to synchotron radiation, which is described by 



dt 



— ^sync-^e 



where the unitless coefficient fogync is given by 



(El) 



(E2) 



ax is the Thomson cross section. Since we consider 
a spherical region extending from 1 kpc to 8 kpc, 
very few of the fermions will be created in the galac- 
tic plane. Therefore, the appropriate magnetic field is 
1 /^G [73], dl], which gives &sync = 6 • 10-43. 

The energy loss of a single fermion due to inverse 
Compton scattering is described by the equation 



where now the unitless coefficient is 



dEe 
dt 



bics — 



iaKNWph 

inrii 



(E3) 



(E4) 



ctkn is the Klein-Nishina cross section, which reduces to 
the Thomson cross section when relativistic corrections 
are negligible. Since this is applicable for scattering with 
CMB photons, 6ics — 5.2 • lO^''^ and is approximately 
independent of energy. (For the parameters with a = 1, 
we have 5.3 • lO"''^ instead.) Since this is two orders of 
magnitude larger than the corresponding value for syn- 
chrotron radiation, we may neglect energy loss due to 
synchrotron radiation. 

Therefore, we calculate the photon energy spectrum 
from inverse Compton scattering with CMB photons. 
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The cosmic microwave background radiation is a black- 
body at TcMB = 2.73 K; therefore the photon density 
per unit energy is 



"-ph(e) 



ph,CMB 



1 



dVde 7r2exp(e/T( 



CMB 



- 1 



(E5) 



where e is the energy of the unscattered photon. For in- 
verse Compton scattering, the number of scattered pho- 
tons per unit energy per unit time produced by an elec- 
tron or positron with Lorentz factor 7 is given by j75| , [76] 

^^(£;,7) = y^ denph(e)aKN(i?,e,7) (E6) 
where ctknI-^', e, 7) is 



(E7) 



and 



G(g, r) = 2qMq) + (1 + 2q){l - q) + 2T^q{l - q) 



4e 



eE 

m 2 



E 



q 



r(me - E) ' 



We have put a subscript on ts to remind us that this vari- 
able measures the time during which the fermion scatters 
against CMB photons. We use the symbol E for the fi- 
nal energy of the scattered photon. The Thomson limit 
corresponds to F <C 1 which is applicable here. By en- 
ergy conservation, only energies E between the following 
values are allowed 



7TOeF 



472 +F 
which we enforce by writing 
d^K 



7meF 



(E8) 



dE dts 



{En) 



de nph(e)crKN(£',e,7) 

e(s,„ax(7,e)-s)e(-B-^„ 



.(7,e)) 
(E9) 



This equation gives the number of photons per unit 
energy per unit time scattered by an electron or positron 
with energy "fm,,. From the fermion energy distribution 
given in equation (16), the corresponding 7 distribution 
is 

F(7) = (ElO) 
V(S2-m2)(l-4m2/m2) 

for 7 between the values 

y(B2 -m2)(l-4m2/m2) 



/max? I mm 



±- 



2mp 



(Ell) 



Averaging d'^N^/dE dt over the 7 distribution gives 



d^N^ 
dE dts 



(E) 



P{l) rfe nph(e)CTKN(£^, e,7) 



.(7,e)). 
(E12) 



Q{E,^U-f,^)~E)eiE~En 




123456789 
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FIG. 7. This plot shows dN^^tot/dE, described by equation 
( E14[ |, evaluated for the first set of parameters (a — 2, ms = 
4 TeV, m„ = 40 MeV). 



This equation gives us the number of photons scat- 
tered per electron (or positron) per unit time; however, 
we require the total number of photons scattered by one 
electron before it loses all of its energy. Properly, we 
should integrate over ts ; this is complicated because 7 is 
a function of ^5. Therefore, we will approximate 



dN^ 



d^N 



dE dE dts 



T, 



(E13) 



where T = l/bicsEe = l/&ics7"ie is the relevant time- 
scale for energy loss. This gives 



dN^ 



1 



— (E) = 

dE hcsme 

■ cTK^iE, e, 7)9(i? 

n 



7m„, 7 JO 

^ij,e)-E)e{E 



de nph(e) 
E, 



(7,e)). 
(E14) 



The equation describes the total number of scattered 
photons of a particular energy, per single electron or 
positron. We have evaluated this equation for the a — 2 
parameters and the result shown in Fig. [7] for energies 
between 1 GeV and 10 GeV. We see that it drops off 
rapidly as a function of energy. 

To find the total number of photons per unit energy 
per unit time, we must multiply by the rate of production 
of the high energy fcrmions, which gives 



dEdt ^ ' ^ 



dt bicsme J^„i„ 7 J a 
aKN{E, e, 7)e(£;„,ax(7, e) " E)e{E - ^n.in(7, e)) 



de nph(e) 



(E15) 



Appendix F: Bremsstrahlung Cross Section 

In this appendix, we derive the cross section for 
bremsstrahlung emission of a cr boson in SS — SS scat- 
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tering. Then 10 relevant tree-level diagrams are shown 
in Fig. [2j Note that the t and u-channel diagrams cancel 
to lowest order in the TOo- — > limit. Therefore, the 
resulting cross section may be smaller than what one 
may naively expect. Let us denote the incoming four- 



momenta as pi and p2, the outgoing momenta of the two 
S particles as and p4, and the outgoing momentum of 
the bremsstrahlung a particle as . The matrix element 
is then 



iM 



^ 

{ml - {p3 -pi)2)(to| 



(Pi+Ps)^) {ml - (p2 -pAY){'m% - (P3 +P5)^) 



Ale, 



A"^ A 



{ml - {p3 

{P3 ^ P4) 



Pi 



)2)(m| 



{P2 - Ps)^) {ml - {p2 - P4)^)(m| 



{Pi - Pb)^) {ml ~ {P2 ~ P4:)'^){ml - (j33 - pi)2) 

(Fl) 



where the last term, in which the momenta p^ and pi 
are switched, represents the contribution of the bottom 
row of diagrams. We will specialize to the center of mass 
frame; we note that the total cross section is a relativistic 
invariant and therefore it is irrelevant what frame it is 
calculated in. Without a loss of generality we write the 
momenta as 



Pi 



ms 



IPil 
2ms 



,0,0,|p. 



(F2) 



P^- ("i5 + |^,0,0,-|p,|), (F3) 



P3^ [ms + 



gL,|P3|sin(03)cos(03), 
|P3|sin(03)sin(03),|P3|cos(03)), (F4) 



Pi = [ms + IP4I sin(6'4) cos{(t)i), 



and 



P5 



IP4I sin(6'4) sin(04), IP4I 008(6*4)) , (F5) 



(v/m2 + |p5|2, Ip^I sin(05), 0, IPsI cos(05)) . (F6) 



Now we turn our attention to the cross section, which 
is given by 



A{Ei+E2)^' 



CTbrcm = / -77^r=^—r=^i'^'^)'^^'^(.Pl+P2~P3-P'l-P5)dLipS 

(F7) 

where the extra 1/2 comes from the two identical parti- 
cles in the final state and dLips is the Lorentz-invariant 
phase space for the final state particles. In particular, 
this is 



dUps = f[ '^'^^ 



i=3 



(F8) 



In the phase space denominators, we may make the 
approximation Ei = E2 = E3 = E4 = ms, and we 
integrate over the 3-momentum delta function, setting 
Pg — — P4— P5. When the S particles are non-relativistic, 
the energy delta function becomes 



ms 



2ms 



2ms 



IPsP (F9) 



Let us call the angle between P4 and pg 6*45. The delta 
function enforces 

\Vl? |P4P IPsP |P4||P5|COS(045) 



ms ms 2ms ms 
which can be solved for |p4| in terms of \p^\ and 6*45 
, , IP51 008(6*45) , 1 



(FIO) 



IP4I 



2 2 
-2\p,\^+A\pj\^ 



(|P5|2COS2(045) 



4ms a/ 



m„ 



1/2 



(Fll) 



We must of course ensure that the result is positive. By 
our choice of coordinates, the d(j)^ integral is trivial; this 
leaves the integrals over 6*4, (p^, 6*5, and d\p^\ to be done 
numerically. This integral is not infrared divergent due 
to the nonzero mass of the a boson. Since the initial 
momentum in the center of momentum frame is pi = 
Vrei/2ms, the above calculation gives a{vrei)- We can 
then average over the relative momentum 



P{Vrel)cr{Vrel) dv 

ret 



(F12) 



using equation (111. 



Finally, we address Sommerfeld factors, which multi- 
ply the cross section and naively can have a large impact 
at low velocities. (Note that this is a multiplicative fac- 
tor in addition to the typical \/v behavior of the cross 
section.) These describe the formation of a quasi-bound 
state during the interaction; the modified cross section is 



CTSo 



■Ka/v 



1 — exp(— Tra/u) 



(F13) 
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For the parameters under consideration, these factors can 
be extremely large, of order lO'^ or 10^. However, it has 
been argued that in this regime the Sommerfeld factor 
given above is unreliable; additional diagrams beyond the 
ladder diagrams implicitly summed in the above equation 
must be taken into account and a proper resummation 



suggests the factors are of order 0(1) to 0(10) [SS]- This 
is supported by some experimental evidence [77], [78], 
including more recent observations at BaBar [73]. As we 
note in the text, even these large Sommerfeld factors (if 
correct) would not be sufficient to produce a detectable 
signal through bremsstrahlung emission. 
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